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Abstract: Aim of this paper is to present a brief and basic introduction 
to polarized beams and the behaviour of a particle's spin in an electromag- 
netic field. It is not intended to give a complete presentation of the spin 
dynamic theory. Instead the most important features are given without 
proof and details not absolutely necessary are left out. Spin depolariz- 
ing mechanisms are discussed and the cures against depolarizing especially 
during acceleration as utilized at COSY are outlined. A list of references 
is given at the end to help the reader to find more deep details on spin 
dynamic theory P^. 



1 Introduction 

The particle spin is an important feature of nuclei and subnuclear parti- 
cles. In general, the particle interactions depend to various degrees on their 
spin states. It is thus not surprising that the experimentalists demand on 
polarized beams at COSY [2] became important rather early. Polarized 
proton beams at COSY are now routinely delivered to internal and ex- 
ternal experiments since 1996. The polarized beams from the ion source 
are pre-accelerated in the cyclotron JULIC [3], injected and accelerated in 
COSY without significant loss of polarization. Imperfection resonances are 
increased in strength by means of vertical orbit distortions, leading to a 
complete polarization reversal (spin flip). Intrinsic resonances are overcome 
by means of a fast ramping air core quadrupole magnet inducing a rapid 
change in tune and therefore preservation of the polarization (tune jump). 

The polarization during acceleration is observed utilizing parts of the 
former EDDA [I] detector system as polarimeter. Only the fast online po- 
larization measurements during the COSY acceleration cycle made possible 
by the EDDA setup allows an efficient correction of the up to 16 first order 
depolarizing resonances encountered. 

Polarized deuterons were accelerated for the first time early in 2003 and 
are also now routinely available for internal and external experiments. Be- 
cause of the lack of depolarizing resonances in the COSY energy range, 



preservation of the polarization is less involved for deuterons as compared 
to protons. However, a polarization measurement technique for vector and 
tensor polarization has to be developed. Also in this case the EDDA detec- 
tor turned out to be the most important and efficient tool to measure the 
polarization even in the accelerating ramp of COSY. 

Polarized proton beams with intensities above 10 10 stored protons with 
a degree of polarization of up to 0.80 are now routinely delivered to internal 
experiments For external experiments polarized protons can be slowly 

extracted from COSY via stochastic extraction without loss of polarization 
by keeping the extraction momentum far enough away from depolarizing res- 
onances and careful adjustment of the betatron tunes during the extraction 
process. 

2 Polarization 

2.1 Polarization of spin-1/2 particle beams 

Along a quantization axis z the spin S of a proton can only be found in 
two eigenstates of S z with eigenvalues of the spin component S z = +K/2 



("up" or "+" state) or S z = -K/2 ("down" or "-" state). The degree of 
polarization is defined as the ratio 



where N± are the number of particles in the two spin states "up" and 
"down". For a 100% polarized beam (P = ±1) the spin states of all particles 
are aligned along a polarization axis. In general three parameters are needed 
to specify the polarization of spin-1/2 particles. Two for the direction of the 
quantization axis and one for the ratio j^- . Thus the polarization of a spin- 
1/2 particle is a vector characterized by a direction and magnitude. This 
type of polarization is therefore called vector polarization. The polarization 
of a bunch of particles is found by an ensemble average. 

2.2 Polarization of spin-1 particle beams 

For a spin-1 system, e.g. deuterons, there are three spin states along a 
quantization axis as depicted in figure ^ 
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Figure 1: Spin states for a spin-1 parti- 
cle. 



Particles having their spin in a m = state have zero vector polarization 
but they are still aligned, i.e. are polarized. Thus the polarization of spin-1 
particles is characterized by vector polarization and an alignment. This is 
called tensor polarization. The vector polarization is given by the ratio 

P = N+ ~ N - (2) 

where N + , Nq, iV_ particles are aligned along the axis with m = 1,0,-1 
respectively. The tensor polarization is determined from 



P, 



N + + iV_ - 2N 
N+ + N + N_' 



(3) 



with -2 < P zz < 1. 

There are several possibilities for vector and tensor polarization. If all 
states were equally populated the polarization would be zero. Some special 
polarization states are: 

• maximum vector polarization P z = ±1 with tensor polarization P zz = 1 

• pure vector polarization (P zz = 0) with P z = ±2/3 

• pure tensor polarization with P zz = — 2 or 1 and P z = 

• P z = ±1/3 with P zz = -1 

The polarized source of COSY is capable to deliver these states as re- 
quired from the experiment. 
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The polarization of a particle P is the quantum average of the spin 
operator S, P = § < ^>\S\ip >. According to the Ehrenfest theorem the 
evolution of the polarization vector (which is called S in the following) is 
governed by a classical equation of motion. 



3 Spin Precession in an Electromagnetic Field 
3.1 Non-relativistic particles 

For non-relativistic non-radiating particles in a magnetic field B the spin 
motion in the particles' rest frame is governed by the equation 

H = ^ B (4) 

that gives the time change of the spin as the torque on the magnetic moment 
jl exerted by the magnetic field. The magnetic moment jl is associated to 
the spin by 

with the rest mass mo and the elementary charge e. The gyromagnetic ratio 
g is 2 for a point-like spin- 1/2 particle. For real particles its deviation from 
2 is expressed by the gyromagnetic anomaly G = {g — 2)/2. For protons this 
value is G = 1.79 and for deuterons G = —0.143. If we insert the magnetic 
moment into eq. |1] the equation of motion of the spin vector is found 

^ = SxQ L (6) 
at 

around the direction of the rotation vector (angular velocity) 

Si = £- a s < 7 » 

which is collinear to the magnetic field, see figure 2. The magnitude of the 
rotation vector gives the time rate of change d<p/dt of the angle 4> covered 
by the spin projetion onto the plane perpendicular to the magnetic field, 



as shown in figured 
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Figure 2: Precession of the spin vector 
around the magnetic held. 



3.2 Relativistic particles 

So far we considered a particle at rest. However in COSY particles are 
moving with a speed close to the velocity of light and also experience an ac- 
celeration when they are bent by the dipole magnets. It turns out that eq. 
is then still valid albeit with a different rotation vector if the electromag- 
netic fields in the laboratory frame are properly Lorentz transformed into 
a co-ordinate system that moves with the particle around the ring (moving 
frame). If for simplicity the particles move in the horizontal plane they expe- 
rience an acceleration that is perpendicular to the velocity and the vertical 
magnetic field (dipole fields of COSY). This is the reason of the Thomas 
precession which contributes an angular velocity in the moving frame. Thus 
the angular velocity £1 in eq.Elwill contain two parts, one component results 
from the Lorentz transformation and the other one from the Thomas pre- 
cession term. The correct form of eq. El for moving particle bears the name 
from Thomas, Bargman, Michel and Telegdi [Thomas — BMT equation) 
and is given by 

^ = 5x n BMT (9) 
with the angular velocity vector 

n BMT = — [(1 + lG)B ± + (1 + G)B\\ - ( 7 G + )/3 x -] (10) 

77lo7 1 + 7 c 

in which the magnetic field B (see figure|3J) in the laboratory frame has been 
decomposed into a componet B\\ parallel to the particle velocity (3 = v/c 
and perpendicular to it, B±. The electric field is E and the kinematic factor 
is given by 7 = ■ „ . 

6 J ' y/T^ 
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Figure 3: Decomposition of the mag- 
netic field into a component parallel B\\ 
to the beam's velocity and orthogonal to 
it B±, in this case a radial and vertical 
component. 



The spin precession has the following properties: 

• An electric field E is equivalent to a magnetic field B=E/c. An electric 
field of 3 • 10 8 V/m has the same strength as a magnetic field of 1 T. In 
general no such electric field is present. In the following electric fields 
are therefore neglected. 

• Comparison of longitudinal (||) and transverse (_L) magnetic fields: 
According to eq. |H] and ^] the spin vector will rotate in longitudinal 
fields by dcj) = VtBMTdt in the time interval dt. Using ds/dt = (5c 
we see that the spin rotates by d<f> = Qbmt(Is/ '((3c) while the particle 
moves the distance ds. Thus the total spin rotation in a longitudinal 
field is 



where we have used that the momentum of the particle is given by 
p = moP'jc. Similarly, for transverse fields the spin rotation is 
determined from 



The spin rotation by longitudinal fields becomes less effective at large 
momenta and the spin rotation by transverse fields becomes nearly 
independent of momentum (7G < 1). 

From eq. ^2 and ^] we can directly calculate that 

• to rotate the spin of 3.3 GeV/c protons around the longitudinal axis 
by an angle tt a longitudinal field with 12.4 Tm is needed at COSY 
and for deuterons: 40.3 Tm (!!) 

• to rotate the proton spin at 3.3 GeV/c in the moving frame by an 



angle tt around the vertical axis one needs a vertical magnetic field 
with / B ± ds = vr^^ = 5.3 Tm, and for deuterons « 118 Tm. 




(11) 




(12) 
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Thus the deuterons spin is about 20 times harder to manipulate with mag- 
netic fields. On the other hand their spin direction is less sensitive to field 
errors. 



4 Spin Precession and Resonances in a Circular 
Particle Accelerator 

We first consider an ideal accelerator or storage ring in which all components 
(dipoles quadrupoles etc.) lie in the horizontal plane. The particles will 
then also move in the horizontal plane on the ideal (circular) path guided 
by the vertical dipole fields. The particle frame rotates with angular velocity 
(cyclotron frequency) 

tic = —B ± . (13) 

m 7 

The spin will precess around the vertical axis and remains constant in the 
straight sections as long as the particle passes the quadrupoles on axis. The 
difference in frequency of spin and velocity precession, eq. El and eq. ITTfl 
yields 

n = vlbmt -n c = jGn c . (14) 

The number of spin precession during one turn in the ring follows from eq. 1141 
and is called spin tune 

v = jG. (15) 

From eq.^Jit is clear that the vertical spin component is preserved while the 
two other components rotate in the moving frame with the angular frequency 
given by eq. El around the vertical axis. As the particle completes one 
revolution, the coordinate system rotates by 2n and the spin has precessed 
around the vertical axis by As a result, in a ring where particles 

experience only vertical magnetic fields the only stable (or invariant) spin 
direction is along the vertical axis. The spin rotation around the vertical 
axis after one turn in the ring can be expressed by the matrix 

cos(2tt7G) sin(27T7G) \ 

10 (16) 

-sin(27T7G) cos(2vr7G) / 

Now consider that at one location in the ring there is a static radial 
magnetic field component (along the x-axis) as shown in figure |3 When the 
particle passes this disturbing magnetic field the spin is rotated around the 
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x-axis by an angle ( given by ea. I12|). The spin rotation around the x-axis 
due to the disturbing field is then given by the matrix 



1 

cos6> 
- sin l 




sin ( 
cos l 



(17) 



The total spin rotation matrix for one turn is then given by the product 
M of ea.UHland eq.EI The stable spin direction outside the disturbing field 
region can be found by determining the eigenvectors of M with eigenvalue 
1, i.e. solving the equation Mu = u. One finds that the vertical spin 
component is: 

y = - Sy (18) 

Vl + tan 2 ^/2)/sin 2 (7r 7 G) 



Clearly, y = S y if there is no disturbing field, i.e. = 0. However polar- 
ization along the vertical axis is completely lost if the spin tune equals an 
integer, v = 7G = k £ N. One can say that the spin motion is in resonance 
with the radial magnetic field: every time the particle passes the field the 
spin vector is bent more and more away from the vertical axis until polar- 
ization is lost. Resonances of this type are called imperfection resonances 
or integer resonances. They arise e.g. if magnets are slightly misaligned or 
if there are vertical orbit distortions. The particles then experience radial 
magnetic fields (e.g. from vertical focussing quadrupoles). Thus it is im- 
possible to operate the accelerator with a polarized beam at energies where 
imperfection resonances occur. However they can be crossed with (almost) 
no polarization loss during acceleration as is discussed below. 

The term "resonance" might become more clear in another view: Each 
time the particle passes the disturbing radial field the spin may point in a 
different direction with respect to the magnetic field if the spin tune is not 
an integer. One can think of a rotation of the radial magnetic field around 
the vertical axis with resonance frequency vrQ,c- While the spin is rotating 
around the vertical axis with spin tune v (not an integer) it rotates with 
frequency 5 = (y — vr)£Ic m the frame of the disturbing field as shown 
in figure 31 If the spin tune is not an integer these rotations are out of 
phase and the disturbing effect of the field averages out. No polarization 
loss occurs. 
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Figure 4: Explanation see text. 



The radial component is according to ea. 1101 given by 



9. 



m 7 



(19) 



where e is proportional to the resonance strength. It is clear that the larger 
the field is the larger the resonance strength will be. Large vertical closed 
orbit errors enhances the imperfection resonance strength. The resonance 
strength can also be interpreted as the width of the resonance since in the 
tune interval (vr — £,vr + e) the angle of f2 with vertical axis is larger then 
45 degrees. On top of the resonance, i.e. 5 = 0, as shown in figureEl the spin 
rotates around the radial direction and on average the vertical component 
is zero and the vertical polarization is lost. 




Figure 5: On resonance the vertical po- 
larization is lost. 



So far we assumed that the particles move, apart from the straight sec- 
tions, on a circle in the horizontal plane. If the particles execute horizontal 
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betatron oscillations the spin rotation is no longer given by the rotation ma- 
trix eq. ED since the path length during one turn of a particle now depends 
on the betatron amplitude and betatron tune Q x . The spin precession is 
modulated by the betatron motion and the spin tune has two additional 
values which are in the case of COSY jG ± (Q x — 2). A horizontal spin 
resonance can occur if 7G ± (Q x — 2) equals an integer. These resonance 
are not important (except when there is coupling between the horizontal 
and vertical betatron motion) at COSY with its vertically polarized beams 
since the horizontal betatron motion is driven by vertical magnetic fields 
of the quadrupoles (see eas. IHland llUj) . In addition to the horizontal beta- 
tron motion the particles also execute a vertical betatron motion with beta- 
tron tune Q y induced by the radial magnetic fields of the vertical focussing 
quadrupoles. In this case eq. El is no longer valid since, besides a possible 
static radial magnetic field, the spin precession experiences, or is driven by a 
force that is evoked by the vertical betatron motion. Similarly, a resonance 
occurs if 7G ± (Q y — 2) equals an integer. These resonances are named 
vertical spin resonances or more frequently, intrinsic resonances. Their 
strength is proportional to the square root of the vertical beam emittance 
(proportional to beam size). Resonances of this type are dangerous in COSY 
because they produce a severe polarization loss if they are not avoided or 
compensated. As in the case of imperfection resonances the spin resonances 
driven by the betatron motion are energy dependent but in addition depend 
also on the vertical tune of the machine. Among these resonances there are 
higher order resonances as well those that can be excited by the synchrotron 
motion of the particles as induced by an rf-cavity (bunched beam). These 
cases play however a minor role at COSY and will therefore not be discussed 
here. There are also resonances that can be excited by longitudinal mag- 
netic fields (solenoids). For this type of resonance the reader is referred to 
the literature. 



5 Resonance Crossing 

Suppose we are away from a single resonance. Then S = [y — z/fj)f2c 7^ 0. 
If the beam is accelerated linearly through the resonance one can write 
v = + a© for the linear variation of the tune with azimuth O. For the 
spin rotation frequency on resonance we have f2j_ = e£lc which shows that 
the change in the spin rotation angle within the "time" AO is A<p = eAO. 
Since the "time" needed to cross the resonance is AO = 2e/a the spin 
rotation angle changes by A(j) = 2e 2 /a during resonance crossing. Thus we 
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can distinguish two regimes: 



• Fast crossing (A</> <C 1) means either a small resonance strength or 
a large crossing speed. There is no time to bent away the spin S from 
the vertical and polarization is conserved. 

• Slow crossing (A<p S> 1) means either a large resonance strength or 
a slow crossing speed. The spin vector rotates much faster around the 
rotation vector Q than the latter moves. The spin vector S follows the 
motion of £1 and is completely reversed. The polarization has changed 
its sign but there is no depolarization (adiabatic spin flip). 

Between these two cases polarization may be completely lost. A quan- 
titative estimate of the final vertical spin component S y compared to the 
initial value before crossing the resonance is given by the Froissart-Stora 
formula 

S v> iniUal = 2e ^ - 1, (20) 

Sy, final 

which includes the two cases of fast crossing, S y ji na i = +S y ^ n i t i a u an d 
adiabatic spin flip, S y j ina i = -S y ^ niUa i. 



6 Polarized Proton and Deuteron Beams at COSY 
6.1 Proton beams 

At COSY vertically polarized proton beams are available in the momentum 
range up to 3650 MeV/c for internal and external experiments. At flat 
top about 1.0 • 10 10 protons with a polarization of more than 80% can be 
achieved. During acceleration to flat top energy up to sixteen depolarizing 
resonances have to be cured for polarized protons. The following graphics 
shows the imperfection and intrinsic resonances in the momentum range of 
COSY. In this figure the imperfection resonances appear as horizontal lines 
since they are independent of tune. The horizontal axis is the fractional 
vertical tune q y (Q y =integer part + q y , for COSY the integer part always 
equals 3). 
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Depolarizing Resonances 




Figure 6: Depolarizing resonances during acceleration of protons in COSY as a 
function of the vertical fractional tunc. Bold: imperfection resonances, semi-bold: 
intrinsic resonances. The shaded zigzag line at q y = 0.62 symbolizes the change in 
working point (not in scale) during acceleration of the beam (tune jump). 

There are six imperfection resonances. A loss of polarization at these 
resonances is avoided by artificially increasing the resonance strength, such 
that the polarization direction is completely reversed (adiabatic spin flip, 
slow crossing) when protons are accelerated across the resonance. A vertical 
correction dipole located at a position with a large vertical beta-function is 
used to increase the orbit distortion that enhances the resonance strength. 
One method to conserve polarization at up to ten intrinsic resonances is a 
rapid vertical tune change (fast crossing). For this a fast tune jump system 
consisting of one air coil quadrupole with a length of 0.6 m and a gradient 
of 0.43 T/m is installed in one arc of COSY. It allows a rapid tune change 
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of about 0.06 within 10 /is. Double crossing of resonances is avoided by a 
slow fall time of 40 ms. Polarization and particle losses due to an emittance 
increase can be kept low during acceleration if the beam position is carefully 
aligned in the acceleration ramp. In addition electron cooling at injection 
energy of COSY has been successfully applied to reduce the beam emittances 
in the transverse as well as in the longitudinal phase space. This increases 
not only the fast crossing (chapter 5) by reducing the resonance strength but 
also lowers beam losses visibile in figure 5 caused by the fast tune jumps. At 
injection the vertical tune is below 3.66. Therefore it is fixed close to 3.62 
during acceleration in order to have enough time for consecutive tune jumps 
(figure EJ). The dynamic tune measurement allows adjusting the tune during 
acceleration as well as the time and amplitude of the tune jumps. Figure [7| 
shows online measurements of the vertical steerer and fast quadrupole cur- 
rents in comparison with the beam current (BCT). The beam is accelerated 
in 2.6 s from injection, 0.294 GeV/c, up to 3.1 GeV/c. Particle losses due to 
the steerer and fast quadrupole actions are less than 10%. Electron cooling 
was not applied in this case. In this example 4 • 10 9 polarized protons are 
accelerated to flat top. 



TeK Run: 100 S/s Hi Res 

[ I -] 




Ch2 1.00V M 1.00s Aux I 1.33 V 

Bil 200mV 

Mathl 20.0mV 1.00 s 



Figure 7: Trace Ml shows the beam current versus time measured with a beam 
current transformer (BCT). Trace 3 represents the current of the vertical steerer 
magnet that excites a total spin flip at the imperfection resonances. The current of 
the fast quadrupole is given by trace 2. 
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Figure |H] shows a polarization measurement versus momentum during 
acceleration to 3.3 GeV/c. The initial polarization is 80% and drops to 75% 
at flat top momentum. The loss of about 6 resonance. 
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Figure 8: Online polarization measurement starting about 1.3 GeV/c. The arrows 
mark the intrinsic resonances which are compensated with almost no polarization 
loss by fast tune jumps. A total spin flip is created with a vertical steerer for the 
imperfection resonances. 

On-line polarization measurements are carried out during acceleration 
with the high precision detector EDDA, designed to measure pp-scattering 
excitation functions during the acceleration of the COSY beam. This de- 
tector is a double layered scintillation hodoscope and in conjunction with a 
CH2 or C fiber target well suited to function as a polarimeter. The great ad- 
vantage of this online measurement is that the behavior of the polarization 
at every resonance is clearly visible in only a few ten minutes, depending on 
the desired measurement statistics. 

6.2 Deuteron beams 

As was pointed out in chapter 3.2 the deuterons spin is about 20 times 
harder to manipulate with magnetic fields. Due to the deuterons G factor 
the lowest imperfection resonance that can occur is jG = 7 corresponding 



14 



to a momentum p = 13 GeV/c which is well above COSYs working range. 
Furthermore there are no intrinsic resonances at COSY in the normal tune 
range 3.5 < Q y < 3.65. Stochastic extraction to the experiment GEM [7j of a 
tensor and vector polarized deuteron beam at 2.4 GeV/c flat top momentum 
has been carried out successfully. 
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